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ABSTRACT 

In 2006 March the Soft Gamma-ray Repeater SGR 1900+14 resumed its bursting activity after about 2 years 
of quiescence. The Swift mission observed the source several times in order to monitor its timing and spectral 
properties. We report on the intense burst "forest" recorded on 2006 March 29 which lasted for ~30 s, when 
Swift was pointing at the source with the narrow field of view instruments. More than 40 bursts were detected 
both by BAT and by XRT, seven of which are rare intermediate flares (IFs): several times 10 42 ergs were 
released. The BAT data were used to carry out time-resolved spectroscopy in the 14-100keV range down 
to 8 ms timescales. BAT and XRT simultaneous data (in the 1-100 keV range) were used to characterize the 
broad-band energy spectra of IFs and verify the results obtained from the BAT-only spectral fits. 

This unique dataset allowed us to test the magnetar model predictions such as the magnetically trapped 
fireball and the twisted magnetosphere over an unprecedented range of fluxes and with large statistics (in terms 
of both photons and IFs). We confirmed that a two blackbody component fits adequately the time-resolved 
and integrated spectra of IFs. However, Comptonization models give comparable good reduced^ 2 . Moreover, 
we found: i) a change of behavior, around ~ 10 41 erg s _1 , above which the softer blackbody shows a sort of 
saturation while the harder one still grows to a few times 10 ergs -1 ; ii) a rather sharp correlation between 
temperature and radii of the blackbodies (R 2 oc kT~ 3 ), which holds for the most luminous parts of the flares 
(approximately for Lt ot > 10 41 ergs -1 ). Within the magnetar model, the majority of these findings can be 
accounted for in terms of thermalised emission from the E-mode and O-mode photospheres. Interestingly, the 
maximum observed luminosity coming from a region of ~ 15 km matches the magnetic Eddington luminosity 
at the same radius, for a surface dipole field of ~ 8 x 10 14 G (virtually equal to the one deduced from the 
spindown of SGR 1900+14). 

We also studied the persistent emission of SGR 1900+ 14 preceding and following the burst "forest" in search 
for transient and permanent timing and spectral variations on short (days/week) timescales. We found that the 
timing/spectral property variations are correlated with the flux and are likely only due to changes of the high 
energy part of the spectrum (above 5 keV), further supporting the twisted "magnetosphere" scenario. 
Subject headings: pulsar: individual (SGR 1900+14) — star: flare — star: neutron — X-rays: burst 



1. INTRODUCTION 

Soft Gamma-ray Repeaters (SGRs) are a small class of 
high energy transient astrophysical sources characterized by 
their emission of short (~ 100ms) bright (10 39 -10 42 erg s _1 ) 
bursts of soft y rays. SGRs have been associated with per- 
sistent X-ray counterparts with l-10keV luminosity of ~ 
10 34 — 10 35 erg s and energy spectra that are generally well 
fit with a power-law model (F ~2). A blackbody compo- 
nent with kT -0.5 keV has been found in SGR 1900+14 
and SGR 1806-20 (Woods et al. 1999; Woods et al. 2004; 
Mereghetti et al. 2005). In contrast to the common short 
bursts, three SGRs have emitted one very powerful long- 
duration burst. These "giant flares" are distinguished by their 
extreme energies (10 44 -10 46 5 erg s _1 ), their hard spectra at the 
onset, and the presence of coherent 5-7 s pulsations at the spin 
period of the neutron star during the decaying tail lasting sev- 
eral minutes. 

Both SGR 1806-20 and SGR 1900+14 were found to spin 



down secularly at a rate of ~ 10 to 10 s s (Kouveliotou 
et al. 1998,1999). Combining SGR periods and their first 
derivatives, and assuming that the B Aip ~ 3.2 x 10 19 (PP) 1/2 G 
relationship generally applied to radio pulsars holds for SGRs 
as well, magnetic field strengths of 10 14 -10 15 G are inferred. 
The latter values provide strong evidence that SGRs are 
highly magnetized neutron stars, i.e., "magnetars" (Duncan 
& Thompson 1992). In the magnetar model, the magnetic 
field is the dominant source of free energy (orders of mag- 
nitudes higher than the rotational energy of the star), power- 
ing the persistent emission through low-level seismic activ- 
ity and heating of the stellar interior (Thompson & Duncan 
1996). When magnetic stresses build up sufficiently to crack 
a patch of the neutron star crust, the resulting "crustquake" 
ejects hot plasma particles into the magnetosphere, which re- 
sults in an SGR burst (Thompson & Duncan 1995). Alterna- 
tively, the short SGR bursts may arise from magnetic recon- 
nection events in the stellar magnetosphere (Lyutikov 2002). 
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TABLE 1 

SGR 1900+14 Swift observation log for pointings carried out during 2006. 
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Obs/Mode 
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Exposure 






(yyyy-mm-dd hh:mm:ss) 


(yyyy-mm-dd hh:mm:ss) 
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BAT/EVENT 


2006-06-10 06:49:02 


2006-06-10 07:09:04 


1202 


00030386001 


XRT/PC 


2006-03-25 22:50:43 


2006-03-27 21:47:47 


49564 


00030386002 


XRT/PC 


2006-03-28 01:05:28 


2006-03-28 23:48:57 


11350 


00030386003 


XRT/WT 


2006-03-29 01:01:00 


2006-03-29 23:57:00 


18592 


00030386004 


XRT/WT 


2006-03-30 01:07:00 


2006-03-30 23:59:58 


20055 
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XRT/WT 


2006-03-31 01:26:23 


2006-03-31 23:59:59 


15931 


00030386006 


XRT/WT 


2006-04-01 01:13:58 


2006-04-03 13:12:31 


45951 
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XRT/WT 
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XRT/PC 
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XRT/PC 


2006-04-12 00:47:49 
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20503 


00030386012 


XRT/PC 


2006-04-13 00:53:42 


2006-04-13 18:41:34 
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00030386013 


XRT/PC 


2006-04-14 22:22:41 


2006-04-14 23:59:56 


1602 


00214277000 


XRT/PC 


2006-06-10 06:54:59 
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"First BAT trigger: 2006-03-25 20: 16:40.28 UT. 
b Second BAT trigger: 2006-03-28 13:50: 15.00 UT. 
c Third BAT trigger: 2006-03-29 01:28:03.99 UT. 
d Fourth BAT trigger: 2006-03-29 02:53:09.46 UT. 
e Fifth BAT trigger: 2006-04-05 06:43:31.36 UT. 
f Sixth BAT trigger: 2006-06-10 06:53:00.8 UT. 

Giant flares likely result from a sudden reconfiguration of 
the star's magnetic field that produces large fractures in the 
crust and propagates outwards through Alfven waves of enor- 
mous power. The ~0.2s initial spike of giant flares marks 
the ejection of the largest part of the total event energy (up 
to ~10 47 ergs) as well as the onset of a relatively large sur- 
face fracture (Schwartz et al. 2005). The detection of deca-, 
hecto-, and kilo-Hertz QPOs in the decaying phases (tails) 
of giant flares have been interpreted in terms of global seis- 
mic oscillations following the fracture formation (Israel et al. 
2005; Strohmayer & Watts 2005). The super-Eddington fluxes 
observed in the tails of these events are thought to be possible 
because of the suppression of the electron scattering cross sec- 
tions in the presence of very strong magnetic fields (Paczyhski 
1992). 

Four confirmed SGRs are known so far, three of which 
are in our Galaxy and one in the Large Magellanic Cloud. 
The Galactic SGRs are located very close to the plane of the 
Galaxy, indicating that SGRs belong to a young stellar popu- 
lation. Furthermore, two of these SGRs are likely associated 
with clusters of very massive stars (Fuchs et al. 1999; Vrba et 
al. 2000; Eikenberry et al. 2001). 

Burst activity in SGRs and, in particular, in SGR 1900+14 
occurs sporadically and is quite diverse in behavior. Burst ac- 
tive phases of SGRs vary both in total energy released and 
duration, tending to be concentrated into relatively narrow in- 
tervals (weeks or months) separated by relatively long peri- 
ods (years) of quiescence. Since its discovery in 1979, when 
bursts where detected three times in three days, SGR 1900+ 14 
remained in quiescence until 1992, when a handful of bursts 
were again detected within a few days (Mazets & Golenetskii 
1981; Kouveliotou et al 1993). In 1998 the source entered an 
unprecedented level of activity during which more than 1000 
burst were recorded within 9 months, and culminated with 
the ~ 400 s long and rather intense 1998 August 27 giant flare 
event (~ 4 x 10 44 ergs _1 ). Burst activity was detected until 



1999. Particularly intersting, are two events detected in May 
and September 1998, during which "bunching" of short and 
long burst were detected (Hurley et al. 1999; Mazets et al. 
1999). Finally, after almost two years of inactivity, on 2001 
April 18, an intense and long (~ 40 s) burst was detected from 
SGR 1900+14 (Guidorzi et al. 2004). The energy released 
by this flare (few x 10 42 ergs) was less than that of the giant 
flare, though larger than that of the most common short SGR 
bursts. Consequently the flare was dubbed an "intermediate 
flare" (Kouveliotou et al. 2001). Indeed, looking back at the 
burst history of SGR 1900+14, it was realized that a hand- 
ful of similar events, characterized by a longer duration (few 
seconds-few tens of seconds) and a higher fluence than those 
characterizing short bursts, were present. These bursts were 
then grouped into the class of intermediate flares likely form- 
ing a continuum in terms of duration and fluences (Olive et al. 
2004). 

A number of broad-band spectroscopic studies of short and 
intermediate flares were carried out in the past by different 
missions. These found that, above 15-20 keV spectra are usu- 
ally well fitted by an optically-thin thermal bremsstrahlung 
(OTTB) with characteristic temperatures ranging from 20 to 
40keV (Aptekar et al. 2001). Laros et al. (1986) and Feni- 
more et al. (1994) studied the cumulative spectral properties 
of 100 bursts from SGR 1806-20 in the 5-200 keV energy 
range finding that the number of the photons in the spec- 
trum of the detected bursts was remarkably stable in spite of 
the large intensity spread (a factor of 50) and that the low- 
energy (below 15 keV) data were inconsistent with the back- 
extrapolation of an OTTB model that provided a good fit to 
the high-energy portion of the spectrum. Qualitatively sim- 
ilar spectral properties were measured during a bright inter- 
mediate flare from SGR 1900+14 using HETE-2 (Olive et al. 
2004) data, which also showed that the OTTB model largely 
overestimated the flux at low energies (<15keV). The broad- 
band spectrum (7-100 keV) was best fitted by the sum of two 
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Fig. 1. — 15-100keV BAT light curves with a time resolution of 1 ms obtained during the burst "forest" of 2006 March 29. 



blackbody (BB) laws. Similar conclusions were reported by 
Strohmayer & Ibrahim (2000) confirming the OTTB limits 
below 15keV by using data from the Rossi X-Ray Timing 
Explorer. 

Feroci et al. (2004) analyzed the 1.5-100 keV BeppoSAX 
spectral properties of 10 short burst from SGR 1900+14, fur- 
ther confirming that the widely used optically-thin thermal- 
bremsstrahlung model provides acceptable spectral fits for en- 
ergies higher than 15 keV but severely overestimates the flux 
at lower energies. A better fit was obtained by means of an 
alternative spectral model such as two blackbodies or a cut- 
off power law. Most of these studies were hampered by one 
or more of the following limitations: relatively poor statistics, 
narrow energy bandpass, poor time resolution, off-axis data, 
and low sensitivity of wide field of view instruments. 

Cumulative analysis of 50 bursts detected by HETE-2 from 
SGR 1900+14 showed the presence of a time lag of about 
2.2±0.4ms between the 30-100 keV and 2-10keV radiation 
bands implying a rapid spectral softening and energy re- 
injection mechanism (Nakagawa et al. 2007). Alternatively, 
the softer spectral component could be reprocessed radiation 
from the harder emission which might be generated near the 
neutron star surface. More generally, the observed time lag 
favors spectral models with (at least) two components. In any 
case the temperatures of the 2BB model does not seem to de- 
pend on either the burst instensity or morphology. 

It has been assessed that burst activity in the SGRs can have 
a persistent effect on the underlying X-ray source (Woods & 



Thompson 2006 and references therein). During the 1998 
burst activity phase of SGR 1900+14, the X-ray counter- 
part became brighter, its energy spectrum was altered, and 
the pulse shape changed dramatically. A similar behavior 
was also recorded from SGR 1806-20 when resumed its burst 
activity phase in 2004 culminating with the 2004 Decem- 
ber 27th giant flare, and later during the 2005-2006 decay- 
ing phase (Woods et al. 2007). It has also been proposed 
that the bursting activity phases are related to the magneto- 
sphere. The correlation between the persistent X-ray flux and 
the spectral hardness observed in several magnetar candidates 
on timescales of years could be accounted for if the evolu- 
tion is regulated by varying "twist" of the magnetosphere, 
as in the case of 1RXS J 170849 -4009 10 (Rea et al. 2005a 
and Israel et al. 2007b) and SGR 1806-20 (Rea et al. 2005b 
and Woods et al. 2007). As discussed by Thompson et al. 
(2002), magnetars may differ from standard radio pulsars be- 
cause their magnetic field is globally twisted inside the star, 
up to a strength of about 10 times the external dipole and, 
occasionally, some of its helicity propagates across the NS 
surface through crustquakes, leading to an impulsive twist up 
of the external field, as well. Such an evolving magnetic field 
is responsible for the different forms of activity and also for 
the glitches. The basic idea of the "twist" scenario is that 
when a static twist is formed, currents flow into the magneto- 
sphere. As the twist angle A0 N s grows, electrons provide an 
increasing optical depth to resonant cyclotron scattering, lead- 
ing to the build-up of a flatter photon power-law component. 
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Fig. 2. — BAT and XRT (WT) light curves obtained simultaneously during 
the burst "forest" of 2006 March 29. Different energy ranges are shown: 
1-4 and -4-10 keV for the XRT (panels XI and X2, respectively), and 15- 
25keV, 25-40 keV, 40-100 keV and >100keV for the BAT (panels Bl, B2, 
B3, and B4, respectively). The XRT light curves were background-subtracted 
and corrected for vignetting, PSF losses, as well as for pile-up effects. It is 
evident from the comparison between XRT and BAT light curves that, on 
average, the IFs are harder than the short bursts, though notable exceptions 
are present. 
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Fig. 3. — Outcoming of the our selection criterion of relying upon a min- 
imum of 4000 photons in each spectrum. For the first three seconds of the 
2006 March 29 "forest" we reported the net counts (upper panel) and the 
integration time (lower panel). 



TABLE 2 

Summary of the average reduced^ 2 and standard deviation for different single- and 
multi-component models obtained for the 729 time-resolved bat spectra and the 8 
bat+xrt spectra of the most intense bursts. 



Spectral 


<xi > 




<xi > 


(T 


Parameters 


Model 


BAT+XRT" 


BAT 




(#) 


Bremss 


4.84 


1.17 


1.71 


1.95 


2 


DiskBB 


2.91 


0.84 


1.17 


0.51 


2 


CompST 


2.41 


0.16 


1.08 


0.42 


3 


CutoffPL 


1.36 


0.07 


1.07 


0.23 


3 


Bremss+BB 


1.33 


0.25 


1.06 


0.28 


4 


CompTT 


0.88 


0.07 


0.99 


0.17 


4 


BB+BB 


0.88 


0.08 


1.01 


0.16 


4 


Bremss" 






0.91 


0.11 


2 



a We used a free scaling factor between XRT and BAT of the order of 10%. 
"Fit carried out in the 15-50 keV range. 

At the same time, returning currents produce an extra heating 
of the star surface and an increased X-ray flux. The few long- 
term X-ray monitoring observations of magnetars collected 
until now are consistent with a scenario in which the twist 
angle steadily increased before the glitch epochs, culminat- 
ing with glitches, sometimes bursts, and a period of increased 
timing noise, and then decreased, leading to a smaller flux and 
a softer spectrum. 

In this paper, we report the results obtained from our anal- 
ysis of the data collected by the Swift Burst Alert Telescope 
(BAT; Barthelmy et al. 2005) and X-Ray Telescope (XRT; 
Burrows et al. 2005) which monitored SGR 1900+14 dur- 
ing a burst active phase in 2006 March (Palmer et al. 2006). 
SGR 1900+14 triggered BAT five times in a few days at the 
end of 2006 March, and once a couple of months later. In 
particular, on March 29 (Romano et al. 2006a) an extremely 
intense emission event (fluence of ~ 2-3 x 10 42 ergs), lasting 
~30s, of short and relatively energetic bursts was recorded 
while both BAT and XRT were simultaneously observing the 
source thanks to an automatic slew performed in response to a 
previous trigger. About 40 single bursts were detected during 
this "forest", or "storm", a handful are long enough (>500ms) 
to be considered "intermediate flares" (IFs). The unprece- 
dented BAT and XRT statistics allowed us to perform both 
average and time-resolved spectral analysis in the 1-100 keV 
band, the first ever on a large sample of intermediate flares 
on ms-timescales. We also studied the persistent emission of 
SGR 1900+14 preceding and following the "forest" in search 
for transient and permanent spectral variations on timescales 
of days and weeks, timescales which have been studied rarely. 
The timing properties of the event are discussed in an accom- 
panying paper (D. Palmer et al. 2008, in preparation). 

2. OBSERVATIONS AND DATA REDUCTION 

Table 1 reports the log of the Swift observations of 
SGR 1900+14 that were used for this work. The BAT data 
of the bursts were accumulated in event mode, with full tem- 
poral and spectral resolution, and archived in sequences la- 
beled by the BAT trigger number (the first six sequences in 
Table 1). The sequence corresponding to the March 29 "for- 
est" is 00203127000. The XRT data of pointed observations 
were archived in the following 13 sequences of Table 1, which 
span from 2006 March 25 to 2006 June 10. Note that the XRT 
data of the "forest" are stored in sequence 0030386003. 

The BAT data were analyzed using the standard BAT anal- 
ysis software distributed within FTOOLS v6.0.5. The ar- 
rival times of events were converted to arrival times in the 
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Fig. 4. — Time-resolved results from the spectral fits obtained in the case of the BB+BB model for the first three seconds of the BAT light curve. The soft BB 
temperature and radius, the hard BB temperature and radius are reported together with the observed flux in the 10-100 keV band and the reduced^ 2 (from the 
top to the bottom panel). 



reference frame of the Solar System Barycenter (SSB) us- 
ing the barycorr task. Mask weighted (i.e. background 
subtracted) light curves were created from event data in the 
15-25, 25-40, 40-100, and 15-100 keV energy ranges at 
1 ms and 4 ms time resolution (see Figure 1 and the four lower 
panels of Figure 2 for the light curves of the March 29 "for- 
est"). Response matrices for BAT spectra were generated with 
the task batdrmgen using the latest spectral redistribution 
matrices. For our spectral fitting (XSPECvl2.2.1) we con- 
sidered the 14-100 keV energy range and applied an energy- 
dependent systematic error vector 1 . The time resolution of 
BAT is 100 microsecs, so pile-up is never an issue under nor- 
mal circumstances (see Barthelmy et al. 2005). No deadtime 
effects have been found for these observations (BAT team; 
private communication). 

The XRT data were first processed by the Swift Data Center 
at NASA/GSFC into Level 1 products (calibrated and quality- 
flagged event lists). Then they were further processed with 
the Heasoft (v6.0.5) to produce the final cleaned event lists. 
In particular, we ran the task xrtpipeline (vO.10.3) apply- 
ing standard filtering and screening criteria, i.e., we cut out 
temporal intervals during which the CCD temperature was 
higher than -47 °C, and we removed hot and flickering pixels 
which are present because the CCD operated at a temperature 
higher than the design temperature of -100 °C due to a fail- 
ure in the active cooling system. An on-board event threshold 
of ~0.2 keV was also applied to the central pixel, which has 
been proved to reduce most of the background due to either 
the bright Earth limb or the CCD dark current (which depends 
on the CCD temperature). 

For our analysis we considered Windowed Timing (WT) 
and Photon Counting (PC) mode data (see Hill et al. 2004 for 

1 http://heasarc.gsfc.nasa.gov/docs/swift/analysis/bat_digest.html 



a full description of read-out modes) and further selected XRT 
event grades 0-12 and 0-2 for WT and PC data, respectively 
(see Burrows et al. 2005). The arrival times of events were 
converted to arrival times in the reference frame of the SSB. 

During sequence 00030386003 (i.e. the "forest") the count 
rate of the outburst was high enough to cause pile-up even in 
the WT mode data. Therefore, to account for this effect, the 
WT data were extracted in a rectangular 40x20-pixel region 
with a 4x20- or 6x20-pixel region excluded from its center 
depending on the peak source intensity, i.e. with a 4x20-pixel 
exclusion region during the first and third snapshot 2 (where 
the observed count rate reached « 350 counts s _1 ), and with a 
6x20-pixel exclusion region during the second snapshot (ob- 
served count rate reaching a 1300 counts s _1 ). The size of the 
exclusion region was determined following the procedure il- 
lustrated in Romano et al. (2006b) and correspond to 30-39% 
(6^1 pixel hole) of the XRT PSF Ancillary response files were 
generated with the task xrtmkarf within FTOOLS, and ac- 
count for different extraction regions and PSF corrections. We 
used the latest spectral redistribution matrices in the Calibra- 
tion Database maintained by HEASARC. 

The PC data show an average count rate of ~ 0.08 counts 
s _1 throughout the entire monitoring campaign, therefore no 
pile-up correction was necessary. We extracted the source 
events in a circle with a radius of 20-pixels (~ 47"). To ac- 
count for the background, we extracted WT events within a 
rectangular box (40x20 pixels), and PC events within an an- 
nular region (radii 85 and 110 pixels) centered on the source 
and far from background sources. 

The energy-resolved light curves during the burst active 
phase is shown in Figure 2 both for the XRT (panels marked 
with XI and X2) and the BAT (panels Bl to B4). The XRT 

2 A snapshot is a continuous pointing at the target. 
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Fig. 5. — Square of radii of the two blackbodies as a function of their temperatures for the whole ~30s BAT dataset. Red filled squares and dark green 
filled circles mark the BB S component for luminosities below and above 3xl0 40 erg s , respectively. Violet filled triangles and blue filled stars are the same 
quantities as above for the BBj, component. Black diamonds marks the Olive et al. (2004) measurements obtained for the 2001 IF. We also plot two representative 
power-laws for R 2 = kT", with a equal -3 and -4 (the latter corresponding to the expected relation for a pure BB.) 



light curves were background-subtracted and corrected for vi- 
gnetting and PSF losses, as well as for pile-up. 

2. 1 . Time-resolved BAT Spectroscopy 

This section refers to the analysis of the BAT dataset 
recorded after the fourth trigger (sequence 00203127000; see 
Table 1). In consideration of the extremely pronounced vari- 
ability of the source during the burst active phase we adopted 
the following strategy for the time-resolved spectroscopic 
analysis: we selected a 4000-count threshold for the accu- 
mulation of each spectrum. This resulted in a set of 729 BAT 
mask weighted (i.e. background subtracted) spectra extracted 
from the event data of sequence 00203127000. The choice 
of the threshold is our best tradoff at optimizing the S/N of 
the spectra, while still closely following the sharpest features 
of the light curve. Figure 3 shows the results of our spec- 
tral accumulation criterion for the first three seconds of the 
BAT light curve. In particular, the choice of 4000 counts per 
spectrum is such that (i) the rapid increase of the first IF is 
well sampled, (ii) the first two short bursts (peaking at about 
3.10 s and 3.35 s) are resolved in time, (iii) 90% of the IF de- 
cay is sampled with a time resolution of <10ms. All spectral 
channels with poor statistics (below 20 counts) were removed 
before fitting. The resulting accumulation time of the BAT 
spectra ranges from 8 ms to about 400 ms. 

Several different spectral models were adopted in the fit of 



the whole sample of 729 spectra in an automatic fashion. For 
each spectrum and model we recorded the values of the best- 
fit parameters, unabsorbed flux, and reduced x 2 - Then, we 
averaged the measured x 1 an d inferred their standard devia- 
tions. The latter quantities can be considered as a qualitative 
indication of the spectral model fit goodness, on average. The 
results are given in Table 2 ordered with the increasing num- 
ber of degrees of freedom in the adopted model. 

Similarly to previous studies reported above and in the lit- 
erature, we tested both single and multi-component models. 
Former models include: an OTTB (Bremss) given that it was 
found to well describe the hard X-ray spectra of SGR bursts, a 
power-law with an exponential cut-off (CutoffPL) which may 
be considered an extension of the simple OTTB, a disk black- 
body (DiskBB) which is a generalization of a thermal spec- 
trum including also geometrical parameters, and different ver- 
sions of a Comptonized spectrum (such as CompST) repre- 
sentative of models involving a thermal spectrum modified 
by mechanisms which preserve the photon number. Among 
double component models we considered a two blackbody 
(BB+BB) able to fit well the broad band energy spectra of 
both short and long SGR burst, an OTTB with a soft black- 
body in order to adjust the OTTB model extrapolation at low 
energy, and a Comptonization model where soft photons are 
up-scattered in a hot plasma taking into account relativistic 
effects (CompTT; Titarchuk 1994) which is, as a first approx- 
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TABLE 3 

XRT WT and BAT spectral fit results of the 7 identified IFs and one short burst (SB) for the BB+BB model. Fluxes are not corrected for the Nu. 



Bin AT H kT R BBs kTh /? BBs F^ 10 ^ F f~ m ™ & d.o.f. 

(10 22 cirT 2 ) (keV) (km) a (keV) (km) a (xl0~ 7 ergs cirT 2 s _1 ) 



IF1 


2.0 ± 0.3 


6.2 ±0.3 


15.9 ±4.2 


12.0 ±0.3 


IF2 


1.9 ±0.3 


5.6 ±0.3 


17.6 ±4.7 


10.4 ±0.3 


SB 


1.2 ±0.3 


5.5 ±0.3 


13.4 ±4.0 


12.1 ±0.4 


IF3 


2.1 ±0.6 


5.4 ±0.3 


17.1 ±5.4 


10.2 ±0.4 


IF4 


2.1 ±0.5 


6.1 ±0.3 


17.4 ±4.9 


11.7 ±0.3 


IF5 


2.8 ± 0.4 


6.1 ±0.3 


17.1 ±4.6 


11.0 ±0.3 


IF6 


2.3 ± 0.4 


5.1 ±0.3 


19.0 ±5.9 


9.8 ±0.3 


IF7 


2.7 ± 0.5 


5.7 ±0.3 


17.1 ±5.1 


11.6 ± 0.3 


a Assumin 


g a distance to the 


source of lOkpc. 








0.1 1 10 

L 4 o(BB soft ) 

Fig. 6. — Time-resolved bolometric luminosity of the BB so f t vs. that of 
the BBhard (in units of 10 40 ). The solid line marks the power-law relation 
L »ft = (L hard )<r with £»=0.70±0.03 (lcr c.l.) obtained by fitting points in the 
0.01-3 L^ {t range. Filled squares represent the measurements by Feroci et 
al. (2004) for short bursts. 

imation, a generalization of both the BB+BB model (in the 
case of saturated Comptonization) and of the CompST one. 

The procedure of increasing new spectral components 
(therefore increasing the number of free parameters) was 
stopped when the F-test probability for the inclusion of an 
additional component became less than 3<x. It is evident from 
Table 2 that small (close to unity) reduced x 2 are obtained, 
for BAT spectra, for one or two component models with free 
parameters > 3. Moreover, small reduced^ 2 are obtained, for 
BAT+XRT joined spectra, only for two component models 
with 4 free parameters. The inclusion of a third component (5 
or 6 free parameters) is not statistically significant. 

The CompTT and BB+BB are by far those that gave the best 
reduced^ 2 (average of 0.99 and 1.01). Moreover, as already 
reported by several authors, while the OTTB model fits the 
spectra well on relatively narrow energy ranges (reduced^ 2 a 
0.9 in the 15-50 keV range), it fails to produce and acceptable 
fit once a broader energy interval is considered (reduced^ 2 » 
1.7 in the 14-100keV energy range). A better description is 
obtained by adding a BB component to the OTTB spectrum 
(reduced^ 2 « 1.07). This additional BB component turned 
out to have a characteristic temperature similar to that of the 



3.5 ± 1.2 


4.68 ±0.15 


103.1 ± 1.8 


0.84 


135 


3.6 ± 1.5 


4.77 ±0.15 


64.4 ± 1.7 


1.15 


153 


2.2 ±0.8 


2.71 ±0.13 


29.1 ±0.9 


0.77 


105 


3.6 ± 1.6 


4.29 ± 0.24 


74.5 ± 1.9 


0.97 


100 


3.7 ± 1.4 


5.49 ± 0.23 


127.8 ±3.0 


0.88 


115 


4.1 ± 1.6 


5.19 ±0.17 


89.8 ±2.4 


1.07 


146 


3.8 ± 1.6 


4.70 ±0.18 


53.5 ± 1.9 


0.73 


117 


3.6 ± 1.2 


4.71 ±0.20 


85.1 ±2.8 


0.77 


116 



softer BB component in the 2BB model. 

The average temperatures and sizes of the BB components 
in the BB+BB during the IFs is <kT soft >= 4.8+0.3 keV with 
<R soft >= 30 + 2 km, and <kT hald >= 9.0 + 0.3 keV with 
<Rhard>= 5.7 + 0.5km, which are in good agreement with 
previous results based on an IF (Olive et al. 2004) and a large 
sample of short bursts (Nakagawaet al. 2007). However, with 
respect to the latter cases, the present Swift data represent the 
best opportunity (so far) to study the above reported spectral 
parameters as a function of time with good statistics coupled 
to fine timing resolution. Figure 4 shows, as an example, the 
results of the automatic spectral fitting in the case of the 2BB 
model for the first 3 s of the high time resolution BAT data, 
while in Figures 5 and 6, we used the inferred values for the 
entire 30 s duration and tried to characterize the main proper- 
ties of the adopted model. 

2.2. Time-resolved BAT+XRT Spectroscopy 

Joint fit of the BAT and simultaneous XRT/WT data of 
the bursting phase associated with the fourth trigger was per- 
formed with a second set of eight BAT and XRT/WT spectra 
(sequences 00203127000 and 00030386003). The eight time 
intervals during which these spectra were accumulated corre- 
spond the occurrence of the seven IFs plus a canonical short 
burst (as a comparison). We report the main spectral parame- 
ters in Table 2 and Table 3. 

The results of the XRT plus BAT fits (see Table 2) con- 
firm those obtained with the BAT in the previous section: the 
CompTT and the BB+BB models gave the best reduced x 2 
(~ 0.9), the OTTB (Bremss Table 2 and hereafter) the worst. 
However, when a BB component is added to the latter com- 
ponent a better solution is obtained (with a number of d.o.f. 
equal to that of the CompTT and the BB+BB models, though 
the reduced^ 2 (~1.3; Table 2) remains too high to be accept- 
able. More generally, the addition of a BB component to the 
other poorly fitted models, such as the CompST and CutoffPL, 
results in a smaller reduced^- 2 , though not yet acceptable and 
with a larger number of free parameters (5). In all cases the 
characteristic temperature of the soft thermal component is 
similar to that of the soft BB in the 2BB model. These find- 
ings, together with similar results obtained by the analysis of 
broad band spectra of SGR bursts from other missions, in 
slightly different (but overlapping) energy intervals, and by 
using different procedures/algorithms, make us confident that 
the Swift fitting results are reliable. 

We note that adopting a different abundance than solar val- 
ues for the absorption multiplicative component (Anders & 
Grevesse 1989; Lodders 2003) does not significantly modify 
the spectral parameters for the joint XRT plus BAT spectral 
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Energy (keV) 



Fig. 7. — The 1-10 keV Swift XRT spectra (PC mode only) of the persistent 
emission of SGR 1900+14 are shown together with the model (BB plus PL) 
residuals. The corresponding Swift data sequences are 00030386001 (black 
lines), 00030386002 (red lines and filled squares), 00030386009 (green 
lines), and 00030386010-6014 (blue lines and filled circles). 




r 



Fig. 8. — The 1, 2, and 3<x confidence regions obtained for the BB temper- 
ature kT and the power-law photon index T for each spectrum during persis- 
tent emission from SGR 1900+14. Note that the temperature kT is consistent 
with being constant, while T is variable, and harder than in quiescence, at 
more than 3<x level. 

fit. 

2.3. XRT Monitoring of the Persistent Emission 

The approximately two-month time span covered by the 
Swift monitoring observations was also used to study the tim- 
ing and spectral properties of the persistent emission. Only 
the PC data were used for spectral analysis in order to rely 
upon a better signal to noise ratio. WT data were also accu- 
mulated and used for timing analysis, after removing bursts, 
in order to rely upon a higher time resolution. The corre- 
sponding Swift data sequences used to carry out the spec- 
tral analysis are 00030386001 (~ 50ks effective exposure 




Norm PL 



Fig. 9. — The 1, 2, and 3cr confidence regions obtained for the BB 
and PL normalizations for each spectrum during persistent emission from 
SGR 1900+14. Similarly to Figure 8 the BB component appears to be ap- 
proximately constant, while the PL normalization is variable and higher than 
in quiescence, at more than 3cr level. 



time), 00030386002 (« llks), 00030386009 (« 29 ks), and 
00030386010-6013 (» 51 ks). In order to minimize the pos- 
sible burst contamination in PC mode, spectra were accumu- 
lated after having excluded ~10s time-interval just before and 
after each good time intervals (GTIs) corresponding to auto- 
matic switches between PC and WT modes. In fact, these 
switches are likely due to high count rate events from cos- 
mic rays and/or from genuine bursts from SGR 1900+14. A 
comparison between the filtered and unfiltered spectra did not 
result in the detection of significant differences, making us 
confident that burst contamination, if any, is negligible. As in 
the previous cases, the data were rebinned with a minimum 
of 20 counts per energy bin to allow x 1 fitting within XSPEC 
(vl2.2.1). We further removed those energy channels with 
poor statistics after background subtraction. 

We assumed an absorbed power-law plus blackbody 
model as inferred with deeper XMM-Newton observations 
(Mereghetti et al. 2006). We left the absorption free to vary 
but forced to be the same in different spectra and obtained 
Nh - (2.3 + 0.1) x 10 22 cm -2 . The results of the fits are re- 
ported in Table 4, while fitted spectra are shown in Figure 7 
together with their residuals expressed in units of cr. There 
is marginal evidence of parameter variability as a function of 
observed flux: the spectrum becomes harder for larger flux 
levels. In order to properly quantify such a spectral variabil- 
ity and its flux dependence, we carried out a bi-dimensional 
fit (contour plot) of the main spectral parameters which has 
also the advantage of better taking into account the possible 
correlations among them. The results of these additional fits 
are shown in Figure 8 and 9. The ellipsoidal shape of the 1, 2, 
and 3<t confidence regions (the innermost region correspond- 
ing to lcr) clearly indicates that parameters (and, more gen- 
erally, the BB and PL components) correlate with each other. 
Additionally, the two-parameter fits confirm the spectral evo- 
lution (c.l. > 3cr) of the persistent spectrum of SGR 1900+14 
as a function of flux: the source displayed harder spectra just 
before (March 28) and after (April 8) the burst "forest" event 
(March 29), and returned back to the pre-event level in mid- 
April. These results are similar to those observed in other 
magnetar candidates, such as SGR 1806-20 among SGRs and 
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1RXS J170849-400910 among AXPs, the main difference 
being in the timescales on which these spectral variations oc- 
curred: years for the latter two sources and days/weeks for 
SGR 1900+14. In the previous cases this trend has been mod- 
elled in terms of "twisting" of the magnetosphere (Thompson 
et al. 2002). 

The WT and PC XRT data were also used to look for the 
-5.2 s pulsations of SGR 1900+14. Given the limited XRT 
statistics in the persistent component and the relatively small 
pulsed fraction of the pulses (~ 15-20% level; semi-amplitude 
of modulation divided by the mean source count rate) we 
used the period value, P=5. 19987+0.00007 s, inferred from 
an XMM-Newton observation carried out on 2006 April 1 
(Mereghetti et al. 2006). We detected the SGR 1900+14 pul- 
sations on 2006 March 29, the same day as burst "forest" 
event, in two time intervals preceding (March 29 from 01:02 
to 01:26 UT) and following (from 04:22 to 04:42 UT) the 
"forest" by approximately 1 hr (in both cases the signal was 
detected at about 3cr c.l.). By merging all the XRT observa- 
tions carried out in March/ April (from 30 March 01:07 to 7 
April 23:20) and by using />=(9.2±0.4)xl0~ n s s" 1 reported 
by Mereghetti et al. (2006), the pulsations could be detected 
at a high significance level. Figure 10 summarizes the results 
of the timing analysis. Both the pulse fractions and shapes in- 
ferred from the first dataset are different from those obtained 
from the second and third datasets, and, more generally, from 
what observed on average since the discovery of this source. 
At the same time, the pulse shapes of the second and third 
datasets are significantly different. The value and shape of 
the the latest dataset is virtually equal (to within the uncer- 
tainties) to those obtained by Mereghetti et al. (2006) about 3 
days after the BAT event. 

2.4. UVOT Observations and Data Reduction 

The Swift UV/Optical Telescope (Roming et al. 2005) is 
a 30 cm modified Ritchey-Chretien telescope co-aligned with 
the XRT. It features a photon-counting micro-channel plate 
intensified CCD detector optimized for UV and optical light 
over the 170-600 nm range; its six broadband colour filters 
sample this range. 

The UVOT was taking data in the V filter (central 
wavelength 5460 A) during the March 29 outbursts of 
SGR 1900+14. No optical counterpart is detected in the re- 
sulting image, accumulated between 2006-03-29 02:39:57 
and 2006-03-29 03:06:00 UT with a total exposure time of 
1539 s. In the following we derived an upper limit to the opti- 
cal emission during the burst "forest". Counts were extracted 
from a 2" radius aperture at the position of the source in the 
XRT; the size of the aperture was chosen so as to minimize 
contamination by a V ~ 18.6 mag star (USNO number 0975- 
14353360) which is only 3" from the XRT position. The 
counts were aperture-corrected up to 6", and used to calcu- 
late a 3cr V-band magnitude upper limit for any counterpart 
of 20.7 mag. Adopting the Akerlof et al. (2000) estimate for 
Ay of ~ 10 mag for the Galactic reddening in the direction 
of SGR 1900+14 we estimate a dereddened upper limit of 
10.7 mag. 

For each UVOT image coincident with an outburst, counts 
were extracted from 4" radius apertures at the position of the 
source in the XRT. The counts were aperture-corrected up to 
6" (optical filters) or 12" (UV filters; the different apertures 
are due to the broader PSF in the UV), and used to calculate 
3cr magnitude upper limits for any counterpart. Using the 
Akerlof et al. (2000) estimate of Ay ~ 10 magnitudes for 
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Fig. 10.— The 1-10 keV Swift XRT lightcurves of SGR 1900+14 from ap- 
proximately one hour before and after the BAT event (2006 March 29; top 
and mid panels), and the first week of April 2006 (bottom panel). In all 
cases we used P=5.19987±0.00007 s, as determined from an XMM-Newton 
observation carried out on 2006 April 1 (Mereghetti et al., 2006). For the 
about one week-long baseline of the April 2006 lightcurve dataset we used 
7 3 =(9.2±0.4)xl0" 11 sr 1 . The pulsed fraction, semi-amplitude of modula- 
tion divided by the mean source count rate, is reported on the right of each 
panel. 

the Galactic reddening in the direction of SGR 1900+14, we 
corrected the upper limits for Galactic extinction assuming 
E(B-V)=3.125. The resulting values are listed in Table 5. 

3. DISCUSSION 

Swift recorded a rather intense and rare series of short bursts 
and intermediate flares (the densest part of which lasted for 
-30 s) from SGR 1900+14 during 2006 March exploiting its 
rapid follow-up pointing capability. This allowed us to carry 
out, for the first time, a detailed study of the average and time- 
resolved spectral properties of this source during such events. 
Moreover, the continuous monitoring of SGR 1900+14 on 
longer timescales, afforded a study of the changes of the spec- 
tral and timing properties of the persistent emission. 

3.1. Burst emission 

We used several different spectral models comprising one or 
two components, to fit both the time-resolved BAT spectra and 
the time-averaged (over a whole single burst or IF) BAT+XRT 
spectra. Interestingly, only two (CompTT and 2BBs) out of 
the five models considered (CompST, CutoffPL, Bremss+BB, 
CompTT and 2BBs) which fit the BAT spectra well, are also 
able to model the BAT+XRT spectra (see Table 2 and Ta- 
ble 3). An optically-thin thermal bremsstrahlung (Bremss) 
does not yield an acceptable reduced^ 2 in either cases (BAT 
and BAT+XRT spectra). This result is consistent with previ- 
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TABLE 4 

xrt pc spectral fit results of the persistent emission of sgr 1900+14 during the four time intervals reported in the text and by assuming the bb+pl 
model. Fluxes are not corrected for the Nu value (which is left to vary but constant among spectra). 



Interval 


kT 


^BBs 


r 


pl-lOkeV 


Y 2 


d.o.f. 


MJD 


(keV) 


(km) a 




(xl0~ 12 ergs crrr 2 s -1 ) 




53819.9519-53821.9082 


0.39 ± 0.02 


3.4 ± 0.5 


2.4 ±0.1 


4.6 ± 0.8 


1.13 


125 


53822.0455-53822.9923 


0.51 ± 0.06 


2.9 ± 0.6 


1.9 ±0.3 


6.3 ± 1.7 


0.91 


39 


53833.0355-53835.9854 


0.48 ± 0.03 


3.3 ± 0.4 


2.1 ± 0.2 


5.0 ± 1.4 


0.96 


86 


53836.0335-53840.0 


0.38 ± 0.02 


4.3 ± 0.8 


2.3 ±0.1 


5.0± 0.7 


0.94 


137 



a Assuming a distance to the source of lOkpc and an absorption hydrogen column of 2.3±0.1 X 10 22 cm 



TABLE 5 

The times of SGR 1900+14 outbursts occurring when UVOT was taking data, with 

THE FILTER IN USE AT THE TIME (CENTRAL WAVELENGTH IN BRACKETS), EXPOSURE TIME OF THE 
IMAGE IN SECONDS, 3tT MAGNITUDE UPPER LIMIT AT THE XRT POSITION. 



Flare time 


Filter 


Exp 


UL 


2006-03-25 23:34:07.748 


UVM2 (2200 A) 


642 


21.3 


2006-03-26 00:59:22.680 


V (5460 A) 


213 


18.9 


2006-03-26 13:30:04.620 


B (4350 A) 


216 


20.8 


2006-03-26 16:50:17.632 


UVW2 (1930 A) 


851 


20.9 


2006-03-28 20:32:31.456 


UVM2 (2200 A) 


261 


20.8 


2006-03-29 02:53:09.46 


V (5460 A) 


1539 


20.7 



ous spectral analyses carried out on broad band spectra of a 
few short bursts and a 4 s-long IF (see, Feroci et al. 2004 and 
Olive et al. 2004). Notably, the Bremss model describes the 
BAT spectra well only when a reduced energy interval, 15- 
50keV, is considered. In the following we discuss in greater 
detail the results inferred from the 2BB model (note that the 
Comp TT model shape is virtually identical to that of the 2BB 
in the limit of saturated Comptonization). Moreover, in con- 
sideration of the difficulties in explaining short bursts and IFs 
in the accretion scenario (they have extreme super-Eddington 
luminosities are in fact observed), we discuss here the impli- 
cations of our results within the magnetar scenario. 

With respect to the former studies of the broad-band spec- 
trum of similar events, the Swift observation of 2006 March 
29 has the advantage of a number of events (>40 between 
short bursts and IFs), therefore largely improving the statis- 
tics over which the detailed analysis can be carried out. The 
impact of the increased statistics and time resolution is ev- 
ident when comparing the Swift results with previous ones. 
On 2001 July 2, the HETE-2 instruments recorded a 3.5 s IF 
from SGR 1900+14 and a time resolved (Af >30ms) spectral 
analysis was carried out in the 2-100keV range (Olive et al. 
2004). They found that the 2BB model fitted the spectra well, 
the higher temperature BB evolving in a manner consistent 
with a shrinking emitting region, the lower temperature BB 
showing a constant radius. They also suggested that, within 
the magnetar model, the 2BB model might be an approxima- 
tion of a more complex multi-temperature spectrum. Their 
Figures 5 and 7 can be easily compared with our Figures 4 
and 5: the Swift data fill in those regions on the kT-R 2 plane, 
which remained unexplored with HETE-2. 

Several new properties can be immediately inferred: a) 
the Swift data populate, almost homogeneously, all tempera- 
tures between ~2 and 12keV, and the HETE-2 measurements 
(black diamonds in Figure 5) can be regarded as a subset of 
them; b) the distribution is such that a sharp edge between the 
populated regions and the rest of the kT-R 2 plane is clearly 
present, with a cut-off (or sharp turn) in the distribution of 



RsBh at kTh ~ 12-13 keV. This sharp edge provides an esti- 
mate of the typical size of the relevant emitting regions (by 
assuming a reference distance of lOkpc): 30-200 km range 
for the BB S component and 3-30 km range for the BB/, one. c) 
there is an additional turn in the BBh component size between 
6 and 16 km, at approximately lOkeV (see Figure 5), and d) 
there is a strong correlation between kT and R BB correspond- 
ing to the brightest phases of the IFs, which is consistent with 
a power-law with index of about -3 (while, were the luminos- 
ity the same for all events, the relationship expected for a BB 
component would be y=x~ 4 , where x - kT and y = R 2 BB ). 

The sharp edge in the distribution of the 2BB parameters 
suggests that a saturation effect occurs. In order to investigate 
this effect further, we divided the data points of the two BB 
components in two samples, below and above a bolometric lu- 
minosity of L = 3xl0 40 ergs _1 (see Figure 5). The luminosity 
threshold was chosen so as to roughly separate the peaks of 
the bursts and flares from the inter-burst/flare time intervals. 
We note that the HETE-2 data lie just above and below the 
flux threshold, but never reach the values corresponding to the 
peaks and tops of the Swift IFs. As we can see from Figure 5, 
as long as each BB component is below ~ 3 x 10 40 ergs _1 , 
their two luminosities are strongly correlated (L{BBh ar d) vs 
L(BB so f t ) is well fitted by a power law with index 0.7+0.3 at 
lcr c.l.). This correlation, which seems indicative of a physi- 
cal link between the two BB components, was already noticed 
by Feroci et al. (2004) who analyzed the time-averaged spec- 
tral properties of 10 short bursts from SGR 1900+ 14 recorded 
by BeppoSAX in April 2001 (filled squares in Figure 6). It is 
apparent that the BeppoSAX measurements are fully consis- 
tent with those obtained with Swift and can be regarded as a 
subsample of the latter. On the other hand, this trend was not 
reported by Olive et al. (2004) for the IF observed by HETE-2 
in 2001 (empty diamonds in Figure 6). Intriguingly, the cor- 
relation no longer holds above ~ 3 x 10 40 erg s _1 for the Swift 
data. Above this value, the luminosity of the BB S component 
reaches a maximum in the (7-14)xl0 40 erg s _1 range, while 
that of the BBh continues to increase up to ~ 3 x 10 erg s _1 . 
Clearly, this previously unknown behavior strengthens the 
idea of a saturation mechanism in the burst emission. 

The emission properties of the higher luminosity flares 
(L > 3 x 10 40 erg s _1 ), combined with the spectral information, 
places some basic constraints on their production mechanism. 
The possible location of the emitting region and the main 
processes likely involved in shaping the spectra have been 
discussed by Thompson & Duncan (1995, hereon TD95). 
Lyubarsky (2002) studied the propagation of a trapped fire- 
ball through the ultramagnetized magnetosphere, and com- 
puted an approximate emission spectrum for it. However, the 
comparison of these predictions with the results from our ob- 
servations or those by Olive et al. (2004) is not straightfor- 
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ward. As already noted by Olive et al. (2004), the "modified 
blackbody" proposed by Lyubarsky (2002) predicts a signifi- 
cant depression below the blackbody of the higher energy tail, 
which is not found in the data (the fit remains good up to en- 
ergies ~ 100 keV). Clearly, more detailed spectral models for 
the burst are required. 

In the following, we will concentrate on new important in- 
formation (mainly concerning the brightest part of IFs) which 
has been inferred from our spectral analysis. In particular, 
we discuss in detail three main findings: i) the existence of 
two blackbody components, the harder one (kT^ ~ 7-1 1 keV) 
having a systematically smaller radius (RsBh ^ 25 km) than 
the softer one (/?bb s ~ 25- 100 km); ii) the existence of a clear 
correlation between the luminosities of the two components 
up to ~ 3 x 10 40 erg s _1 , above which value the luminosity of 
softer blackbody shows signs of saturation, while that of the 
harder blackbody still grows up to a few times 10 41 erg s _1 ; iii) 
the existence of a correlation between temperature and black- 
body radius, which holds for the most luminous parts of the 
flares (approximately for L tot > 4 x 10 40 erg s _1 ). 

The first finding shows that different thermal components 
originate from different regions around the NS, ranging from 
its surface (at R < R^s) to well up in the magnetosphere, at 
a height of several stellar radii. The maximum observed flare 
luminosity is ~3xl0 41 erg s _1 , attained by the hard blackbody 
component at an effective temperature of ~10keV and radius 
of ~15km. Interestingly, this matches well the magnetic Ed- 
dington luminosity, Z-Edd.B. at that same radius, for a surface 
dipole field ^8xl0 14 G (this value of the B-field is very close 
to that deduced from the spindown rate of SGR 1900+14), 
where (cfr. Paczynsky 1992, TD95): 
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If radiation originated from a trapped hot fireball this would 
be in line with the discussion of Thompson & Duncan (1995, 
2001), according to which the radiative efficiency of a mag- 
netic confined fireball never exceeds (to within a small fac- 
tor) the magnetic Eddington flux. In fact, ablation of mat- 
ter from the NS surface would increase the scattering depth 
in a super-Eddington radiation field, thus providing a self- 
regulating mechanism. The radius and temperature that we 
inferred at maximum burst luminosity are also in good agree- 
ment with the prediction for the emission coming from a 
trapped fireball (TD95, cfr. their discussion for the best fit 
temperature and radius in the case of SGR 1806-20). 

Most remarkably, we find that the luminosity of the soft 
component can be larger than 10 40 ergs~' out to R~100km. 
For R>30-40km this value exceeds the magnetic Eddington 
luminosity at the corresponding radius (for a dipole field ge- 
ometry) and suggests that magnetospheric confinement plays 
an active role out to those distances, with magnetic stresses 
balancing radiation forces. If this were the case, the satura- 
tion of the soft component at 10 41 erg s _1 would be related 
to the magnetic field strength at the relevant radii and its max- 
imum ability to retain the trapped fireball matter which is sub- 
ject to very high radiation pressure. A detailed investigation 
of the possible equilibrium configurations is beyond the scope 
of this paper and will be addressed in a future work. 

Our second finding [point ii)] implies that, in the luminosity 
range 10 40 — 10 41 ergs -1 , the total radiation energy is divided 
almost equally between the two components. Combining the 
information of Figures 5 and 6 we can see that, for a given 



luminosity < 10 41 erg s _1 , there are two separate emission re- 
gions: a smaller and hotter one, whose radiating area suggests 
emission from (part of) the NS surface with a relatively high 
effective temperature (10-12keV), plus a second, possibly 
magnetospheric region, with significantly larger emitting area 
and lower effective temperature (3-7 keV). The hotter compo- 
nent attains the highest luminosities, i.e. > 10 41 erg s _1 , as its 
radius grows slightly and the effective temperature decreases. 
On the other hand, the colder component is characterized by 
a minimum blackbody radius of ~ 25 km and maximum tem- 
perature of ~ 7 keV. 

A possible interpretation involves the different way in 
which photons with ordinary (O) or extraordinary (E) polar- 
ization mode (a properties introduced by the presence of a 
birefringent medium, such as the magnetic field; Meszaros 
et al. 1980) propagate across the magnetosphere. Since the 
scattering cross section of E-mode photons is much reduced 
in presence of strong magnetic fields, E-mode photons have 
a scattering photosphere which is located much closer to the 
NS than that of O-mode photons. On the other hand, in su- 
percritical magnetic field (B > Z? QED ), E-mode photons have 
a non-negligible probability of splitting (and switching to the 
O-mode). The probability is a strong function of energy, so 
that E-mode photons of high energy cannot travel far from 
their scattering photosphere before splitting into O-mode pho- 
tons. On the other hand, O-mode photons will be entrained 
with streaming electrons and baryons and advected out to 
where their scattering optical depth becomes r( O) ~ 1 (TD95; 
Lyubarsky 2002). As a consequence of Compton scattering 
(and, with lower efficiency, photon merging), a fraction of O- 
mode photons can also switch back to the E-mode. There- 
fore, the two modes are effectively coupled and are advected 
at comparable rates as long as mode switching is efficient: an 
approximately equal distribution of energy in the two modes 
is thus to be expected (TD95; Lyubarsky 2002). This may 
suggest that the two observed spectral components reflect the 
population of photons in the two polarization modes and, thus, 
the regions from which they are emitted. 

We note that the minimum radius of the cold blackbody 
almost corresponds to the maximum radius of the hot com- 
ponent, the distribution of points in the R 2 vs. kT plane be- 
ing continuous across this point (see Figure 5). This sug- 
gests the existence of a narrow zone of separation between 
the two emission regions. In order to further characterize 
this zone, we studied the R 2 vs. kT distribution for each of 
the seven longest IFs and for one short burst we already dis- 
cussed in previous sections (see also Table 3 and Figure 11). 
Remarkably, the transition between the two BB components is 
more evident for the brightest flares (IF1 and IF7 in Figure 11) 
where the upper left edge of the parameter distribution of the 
hard component is confined to smaller radii (R < 20 km) and 
correspondingly higher temperatures (kT > 8 keV). The soft 
component seems less dependent on the flux than the hard 
one, its minimum size being of the order of > 25-35 km in 
all IFs. Qualitatively, this behavior is in agreement with the 
above scenario. In fact, a larger brightness can be reached 
when the burst energy is trapped closer to the NS surface, 
where the stronger magnetic field allows a larger flux (of E- 
mode photons) to be released. Moreover, the sharp and nar- 
row gap between the emitting regions of the two BBs, at a 
radius of (~ 20-30 km), can be interpreted as the signature 
of the presence of an E-mode splitting photosphere (TD95); 
this represents the surface (at radius R = R QED ) below which 
the magnetic field is supercritical (B > B QED ), photon splitting 
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is efficient and the two polarization modes can remain cou- 
pled. E-mode photons can stream freely to the observer from 
the splitting photosphere (at radius R QED ~ 20-30 km), which 
then would naturally define a maximum size for the emission 
region of radiation in this polarization mode. On the other 
hand, O-mode photons could not originate from regions in- 
ternal to the splitting photosphere, as their scattering optical 
depth is still much larger than unity within that volume. 

If this scenario is correct, it implies that the emission region 
of E-mode photons can range from the base of the trapped 
fireball (near the NS surface) to a slightly higher region in the 
magnetosphere, but well within the splitting photosphere. In 
this case, part of the internal energy of the coupled E and O- 
mode photons is lost in the E-mode photosphere (adiabatic) 
expansion from near the NS to the splitting photosphere. This 
degrades the original spectrum towards lower temperatures 
and larger radii, up to R QED . On the other hand, O-mode pho- 
tons might either be be released near the splitting photosphere 
or advected to much larger radii, depending on the value of 
their optical depth at that location (see TD95). In any case an 
improved treatment of the above effects would be required in 
order to develop a self-consistent model. 

Finally, our third finding is the apparent correlation between 
the surface and temperature of the spectral components at dif- 
ferent luminosities. This behavior does not have a clear inter- 
pretation. As shown in Figure 5, the top part of the data strip 
is well fit by a power-law S oc r~ B 3 , where S is the radiating 
blackbody surface. Since T.T* B = L, this implies r BB oc L. 
Therefore, along lines of constant luminosity, the radiating 
surface would scale like T~ 4 , i.e. steeper than the data points 
(on the edge of the allowed region in the kT vs R 2 plane). The 
Zr 3 =const relation that implies a roughly constant number 
of emitted photons per unit time. 

At the highest luminosities (> 10 41 ergs~') of the flares, 
where the softer blackbody saturates, the observed correlation 
involves both components. It is tempting to attribute this to 
the fact that, in this luminosity range, Comptonization, and/or 
adiabatic losses play a dominant role in the formation of the 
emerging spectrum, both processes being characterized by 
conservation of the photon number (photon splitting should 
have only a minor impact). This is in agreement with the 
prediction by TD95, that at effective temperatures less than 
~ lOkeV photon splitting would not be efficient in maintain- 
ing a pure Planck spectrum for both modes. However, if the 
photon chemical potential is sufficiently small compared to 
the photon temperature, deviations from a Planck spectrum 
are expected only at low energies (TD95; Lyubarski 2002), 
which we do not observe given they are masked by absorp- 
tion. 

On the other hand, at luminosities < 10 41 ergs~', the 
two blackbodies approximately share the same luminosity, 
therefore the colder one is emits a larger number of pho- 
tons. This suggests that efficient photon number-changing 
processes (such as photon splitting or double Compton scat- 
tering) may play a key role in the formation of the spectrum. 
Indeed, we note that the photon number ratio is roughly the 
same as the temperature ratio in this case, since the luminosi- 
ties are approximately equal. As this ratio is peaked around 
2 and extending from 1 .5 and 4, it appears to be qualitatively 
compatible with the combined effects of photon splitting and 
double Compton scattering (cfr. Lyubarsky 2002). 

Finally, we briefly comment on the sharp turn in the hard 
component distribution at around kT~10keV and subsequent 



cut-off at ~12+14keV which, together with the ET -const re- 
lation, are among the most distinctive properties of Figure 5. 
Among others, these features clearly imply that no thermal 
component with characteristic temperature above the cut-off 
is detected, though the BAT energy interval extend up to at 
least lOOkeV. Within the discussed scenario it is tempting 
to account for the above properties in terms of photon split- 
ting process and, in particular, its efficiency as a function of 
(mainly) flux and energy. In fact, the turn at kT~ 1 keV marks 
the drift of the parameter distribution from the ZT 3 =const 
relation, identifying the passage between photon number- 
conservation (such as Comptonization) and photon number- 
changing (photon splitting) emitting processes. As reported 
by TD95, the minimum temperature at which photon splitting 
is efficient in maintaining a nearly Planckian spectrum for E- 
mode photons, is ^10keV. 

3.2. The persistent emission 

The Swift XRT monitoring data allowed us to study the evo- 
lution of the persistent emission of the SGR 1900+14 on rel- 
atively short timescales starting from one day before up to 
one an half month after the BAT event. We clearly detected 
a spectral evolution in the higher energy end of the l-10keV 
spectrum as a function of flux, the PL component becomes 
flatter for larger fluxes. The changes in the persistent com- 
ponent are temporary and last for approximately a couple of 
weeks, while we detected variability on daily timescales. It is 
also interesting to note that the spectral results we obtained on 
2006 April 8 (see P3 interval in Table 4 and the green lines in 
Figures 7, 8 and 9) are consistent, to within the uncertainties, 
with those reported by Mereghetti et al. (2006) from a deeper 
ToO XMM-Newton observation carried out on 2006 March 1 . 

Both AXPs and SGRs have shown a correlation between the 
X-ray flux and the spectral hardness (see, for example, Woods 
et al 2007; Rea et al. 2005a); this has been explained in terms 
of the onset of a "twist" in the magnetosphere. As discussed 
by Thompson et al. 2002, magnetars may differ from standard 
radio pulsar because their magnetic field is globally twisted 
inside the star, up to a strength of about 10 times the external 
dipole and, occasionally, this field is expected to fracture the 
crust and lead to a twist of the external field. The basic idea is 
that when such a twist is formed, currents flow into the mag- 
netosphere. As the twist angle A0 N s grows, electrons provide 
an increasing optical depth to resonant cyclotron scattering, 
leading to the build up of a flatter photon power-law com- 
ponent. At the same time, larger returning currents produce 
an extra heating of the star surface and increased X-ray flux. 
Therefore, both increased activity and glitches are expected 
to be associated with erratic fracturing of the crust. Obser- 
vations obtained until 2003 were consistent with a scenario in 
which the twist angle was steadily increasing before the glitch 
epochs, culminating in glitch(es) and period of increased tim- 
ing noise, and then decreasing, leading to a smaller flux and 
a softer spectrum. The decay time of the global twist is given 
by (Thompson, Lyutikov & Kulkarni 2002) 

(2) 

where A0 N _ S is the net twist between the north and south mag- 
netic hemispheres, B po i e is the magnetic field strength at the 
pole(s). By using the inferred persistent unabsorbed luminos- 
ity of SGR 1900+14 (~10 35 erg s^ 1 ), the observed time inter- 
val, of the order of < 1 month, over which the spectral param- 
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eters change and return to the pre-burst "forest" event, and 
assuming a standard size for the neutron star surface and a 
10 14 G < flp ie G < 10 15 we infer a twist angle in the few 10 3 
to few 10~ 2 radians interval, which is similar to that inferred 
for the bursting AXP IE 2259+586 (~ 10~ 2 rad; Woods et al. 
2004). 

We regard as unlikely that the hardening and softening of 
the power-law component (in the l-10keV band) is due to 
a change of the flux and/or photon index of the additional 
power-law component detected by INTEGRAL in the 20- 
200 keV band flux of several AXPs and SGRs (Kuiper et 
al. 2006; the SGR 1900+14 emission above 20keV was re- 
discovered also by PDS on board BeppoSAX; Gotz et al. 2006; 
Esposito et al. 2007). This is implied by the finding that the 
spectrum of SGR 1900+14 can be accounted for by only one 
power-law component covering both the soft and hard X-rays 
(from~l to -200 keV) 

We monitored the pulse shape and pulse fraction of the 
~5.2 s spin period modulation through a timing analysis of the 
persistent component. Large changes were detected as a func- 
tion of flux (see Figure 10): pulses appear to be less pulsed but 
more structured (to within the uncertainties) and with smaller 
relative amplitude for increasing fluxes. Such a flux depen- 
dence of the pulse properties was already detected in other 
SGRs and AXPs and related to changes in the persistent emis- 
sion component. Among SGRs, SGR 1806-20 showed a de- 
crease in the pulsed fraction together with the appearance of 
additional peaks in the pulse profile after the giant flare of the 
2004 December 27 (Rea et al. 2005b, Woods et al. 2007). 
More recently, a similar behavior was observed in the tran- 
sient AXP CXOU J 1647 10.2-4552 16 in correspondence to 
the onset of an outburst (Israel et al. 2007a). Given the almost 
constant properties of the BB component, variations are likely 
due to the PL component which is supposed to originate into 
the magnetosphere. In particular, The pulse shape changes 
may indicate an evolution of the geometry of the magneto- 
sphere. Alternatively, the magnetosphere may randomize the 
outcoming pulsed photons, temporarily decreasing the pulsed 
fraction. 

4. CONCLUSIONS 

On 2006 March 29, after several days of burst activity, the 
Soft y-ray Repeater SGR 1900+14 displayed a series of short 
bursts and intermediate flares lasting for about 30 s, during 
which all the Swift instruments were pointed to the source. 
The good statistics, the fine time resolution (4 ms) and the 
wide energy coverage (0. 2-100 keV) allowed us to carry out 
a detailed study of the timing and spectral properties dur- 
ing both the bursts and persistent emission. Large varia- 
tions were detected and their evolution studied over a range 
of timescales. Although not all the observational properties 
can be accounted for by the current magnetar scenario, the 
large majority of them are in overall good agreement with it. 
In particular, we found: 

• a break, around ~ 10 41 ergs _1 , in the known corre- 
lation between the luminosities of the two blackbody 
components which fit well the BAT time-resolved and 
BAT+XRT integrated spectra. Above this value the 
softer blackbody shows signs of saturation, while the 
luminosity of the harder blackbody still increases up to 
~ 3 x lO^ergs- 1 ; 

• the existence of a correlation between temperature and 
emitting surface of the blackbodies (R 2 oc kT~ 3 ), which 



holds through the most luminous phases of the flares 

(Lt ot > 104'ergs- 1 ). 

The above two findings together can be interpreted in terms 
of the two populations of photons from the O- and E-mode po- 
larizations modes which are foreseen in the magnetar scenario 
and the regions from which they escape: larger and colder the 
O-mode photon region, harder and smaller the E-mode one. 

Moreover, we note that the maximum observed luminos- 
ity from flares (~ 3 x 10 41 ergs -1 ) is attained by the harder 
blackbody component at an effective temperature of ~10keV 
and radius of ~15km, which is similar to the inferred mag- 
netic Eddington luminosity at the same radius, for a surface 
dipole field ^ 8 x 10 14 G (as deduced from the spindown of 
SGR 1900+14). 

Finally, based on our analysis, we do not see any signifi- 
cant difference between IFs and short bursts. Both IFs and 
short bursts seem to form a continuum in terms of spectral 
properties (besides duration and fluence). 

As far as the persistent emission properties are concerned, 
we found: 

• a clear correlation between the PL photon index F and 
the persistent flux, such that flatter PL components cor- 
respond to larger fluxes. The source returned to the 
usual spectral shape after approximately one week from 
the burst "forest" event; 

• changes in the modulation of the 5.2 s spin period, with 
the pulse profiles becoming more structured (to within 
the uncertainties) for increasing persistent flux. The 
pulse properties returned to their usual state within one 
week, as well. 

The above results are not dissimilar from those observed 
in other magnetar candidates, although this is the first time 
that variations are detected on timescales of an hour (for the 
pulses) and days (for the spectra). These findings are qualita- 
tively in agreement with the expected properties of a "twisted" 
magnetosphere. 

More generally, the unprecedented high quality of the Swift 
data allowed us to gather new important insights into the emis- 
sion mechanism properties, at a detail level which had not 
been attained before. We hope that findings like the ones pre- 
sented in this paper will stimulate a more detailed theoretical 
treatment of the physics of SGR and AXP bursts. 
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